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Abstract. The statistical properties of circumstellar disks around young stars are important for 
constraining theoretical models for the formation and early evolution of planetary systems. In 
this brief review, I survey the literature related to ground-based and Spitzer-based infrared (IR) 
studies of young stellar clusters, with particular emphasis on tracing the evolution of primordial 
("protoplanetary") disks through spectroscopic and photometric diagnostics. The available data 
demonstrate that the fraction of young stars with optically thick primordial disks and/or those which 
show spectroscopic evidence for accretion appears to approximately follow an exponential decay 
with characteristic time ~2.5 Myr (half-life ~ 1.7 Myr). Large IR surveys of ~2-5 Myr-old stellar 
samples show that there is real cluster-by-cluster scatter in the observed disk fractions as a function 
of age. Recent Spitzer surveys have found convincing evidence that disk evolution varies by stellar 
mass and environment (binarity, proximity to massive stars, and cluster density). Perhaps most 
significantly for understanding the planeticity of stars, the disk fraction decay timescale appears 
to vary by stellar mass, ranging from ~1 Myr for >1.3 M stars to ^3 Myr for <0.08 M brown 
dwarfs. The exponential decay function may provide a useful empirical formalism for estimating 
very rough ages for YSO populations and for modeling the effects of disk-locking on the angular 
momentum of young stars. 

Keywords: stars: pre-main sequence, planetary systems: protoplanetary disks, planetary systems: 
formation 
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INTRODUCTION 

For the 2nd Subaru International Conference on "Exoplanets and Disks: Their Forma- 
tion and Diversity", the goal of my introductory talk was to succinctly summarize some 
recent observational constraints on the parameters describing circumstellar disks rele- 
vant to planet formation. Rather than give short shrift to many subtopics within too few 
pages, I have elected to focus this manuscript on only one aspect of disk evolution: pri- 
mordial disk lifetimes. For reviews on recent results on other aspects of circumstellar 
disk evolution, I refer the reader to the other reviews in this volume that were presented 
at the Kona meeting: especially those on disk geometry and coronagraphic polarime- 
try by M. Perrin, observations of debris disks by A. Mora-Martin, dynamical theory of 
planets, planetesimals, and dust by M. Wyatt, imaging debris disks by P. Kalas, mod- 
eling circumstellar disks by S. Wolf, spectroscopy of disks by M. Goto, characterizing 
gas and dust around main sequence stars by C. Chen, and dust growth in disks by H. 
Tanaka. For more exhaustive recent reviews on circumstellar disks around young stars, 
I refer the reader to the book by Hartmann [Q]], and the reviews by Beckwith et al. 10], 
Hillenbrand [H], Hillenbrand [4], Meyer et alfll], Monin et al. [@], and Alexander n[]T 



Investigations of circumstellar disks and extrasolar planets have experienced explo- 
sive (and seemingly parallel) growth over the past two decades. Our knowledge of ex- 
trasolar planets has been driven mostly by enhanced techniques of optical spectroscopy 
(e.g. Doppler surveys), and more recently, photometry (e.g. transits, microlensing). Our 
knowledge regarding circumstellar disks orbiting young stars has been driven mostly by 
infrared and millimeter photometry and imaging, and more recently, spectroscopy (with 
obvious contributions to understanding accretion from optical spectroscopy). Protoplan- 
etary disks set the initial conditions of planet formation, and to understand the diversity 
of planetary systems, we need to understand the physics, chemistry, and evolution of 
primordial disks. 



PRIMORDIAL DISKS 

The existence of primordial disks around young stars was originally inferred through 
the spectroscopic evidence for accretion among classical T Tauri stars (<2 M pre- 
main sequence stars) and evidence of circumstellar dust structures, later confirmed to be 
geometrically distributed as disks [see review by and references therein]. Classical 
T Tauri stars have typical isochronal ages of <5 Myr, typical accretion rates of ~10~ 7 - 
1(T 9 Moyr- 1 H], disk masses of ~5 x 10 3 M© (~0.5 dex dispersion), and median 
disk/star mass ratios of ~0.5% [9]. Among the classical T Tauri star populations in star- 
forming regions are weak-line T Tauri stars - pre-MS stars lacking spectroscopic of 
accretion and/or photometric evidence for circumstellar dust. The disk masses inferred 
for weak-lined T Tauri stars from submillimeter observations are almost always < 10~ 3,5 
M [|9j]. Integration of the minimum-mass solar nebula model over the range of orbital 
radii for planets in our solar system leads to a disk mass of ~10~ 2 M , [i.e. ~2x the 



typical disk mass inferred for classical T Tauri stars;|10D, indicating that by the time low- 
mass stars evolve to the weak-line T Tauri phase their disk surface density is likely to be 
at least an order of magnitude too low to form planets like Jupiter and Saturn. By ages 
of ~10 Myr, samples of typical solar-type stars appear to have less than <tens of M,Earth 
of gas within a few AU, and <few M £art /, of circumstellar gas at orbital radii of ~ 10-40 



AU fl 1 IN . Early in the protostar's life, the disk (and protostar) are fed by infall from a 
vast molecular envelope associated with cloud cores. The molecular clouds forming em- 
bedded star clusters appear to disperse within ~3 Myr of star formation II 1211 . effectively 
removing the source of dense gas feeding primordial circumstellar disks. The mech- 
anisms for removing reservoirs of circumstellar gas are numerous, including viscous 
accretion radially inward toward the star (with some material subsequently ejected out- 
ward via jet), outward viscous decretion (through conservation of angular momentum), 
accretion into planets, photoevaporation by the central star, or even photoevaporation 
by a neighboring star Q_l], LZO - Given the numerous mechanisms for dispersing primordial 
disks, it is perhaps unsurprising that we find that typical disk lifetimes are similar to the 
brief timescale when the protostar is immersed in a sea of dense molecular gas. 

While the longevity of the Sun's protoplanetary disk is unknown, there are weak 
constraints. Detailed modeling of the geophysical, thermal, and rotational evolution of 
Saturn's outermost large, regular satellite Iapetus by Castillo-Rogez et al. [13] requires 



that Iapetus accreted most of its mass within ~3.4-5.4 MyiQ after the formation of the 
solar system [defined by the formation of Ca-Al-rich inclusions 4,567.2 ± 0.6 Myr ago; 
[l5ll . This suggests that Saturn itself had accreted the majority of its mass [which is 
at least ~67-80% hydrogen & helium; LL6J] from the Sun's protoplanetary disk within 
~3.4-5.4 Myr. Variations of the 'Nice model' by Thommes et al. M 1 VL 1 1 811 can form 
the gas giants Jupiter and Saturn through core accretion, form Uranus and Neptune as 
'failed' gas giants in the ~5-10 AU zone, and scatter them to near their current orbital 
radii - all within 5 Myr - for models where the disk surface density is roughly an order 
of magnitude higher than the canonical Hayashi lfl9ll minimum-mass solar nebula. These 
results suggest that the formation of giant planets is at least plausible within the ~10 6 ~ 7 
yr lifetime of gas-rich protoplanetary disks, given the observed physical properties of 
observed primordial disks. 



PRIMORDIAL DISK LIFETIMES 

Photometric and spectroscopic observations appear to be telling us that the circumstellar 
disks of young stars undergo radically divergent evolutionary paths at a very young 
age. Early ~3-4 /im surveys suggested 1I20I1 . and recent Spitzer surveys of nearby young 



stellar groups have confirmed 112 U 12211 . that primordial accretion disks around lower- 
mass stars tend to last longer than around higher- mass stars. However within a given 
mass range, there appears to be quite a dispersion in disk lifetimes. Shorter disk lifetimes 



have been demonstrated for stars which are higher in mass [e.g. |22J], in multiple systems 



fe.g. I23L 12411 . an d which are in the immediate vicinity (<0.5 pc) of O-type stars [e.g 



251 126H . or are in denser stellar clusters [e.g. 27]. Hence there appear to be multiple 
variables affecting primordial disk lifetimes. Although the well-cited disk survey of 
Haisch et al. [28] pointed towards a maximum disk lifetime of 6 Myr, we now know 
of many convincing examples of older accretors: multiple examples in the ~6-10 Myr- 
old 77 Cha, TW Hya, and 25 Ori groups, the accretor PDS 66 in Lower Cen Crux (stellar 
age ~7-17 Myr), and the unusual binary StHa 34 (~8-25 Myr) 0]. 

A modern version of the "Haisch-Lada 2 plot" (disk fraction vs. sample age) is shown 
in Fig. [TJ While the observational definitions of what constitutes a star with a probable 
primordial accretion disk can vary slightly from study to study, this plot represents a best 
effort to summarize the situation with the data available. As demonstrated by Haisch 



et al. [29] from the sample of known T Tauri stars in Taurus [|43|], 3-4 /im excesses 
trace spectroscopically identified classical T Tauri stars ~100% of time, while ~2 /im 
excesses trace classical T Tauri stars only ~70% of the time. Similarly, Silverstone et al. 
[14441 detected 3.6 /im excess emission solely around previously known classical T Tauri 
stars among a large sample of FGK-type stars with ages ~3 Myr to ~3 Gyr. With these 
findings in mind, I consider stars to be accretors whether they are (1) spectroscopically 
identified as classical T Tauri stars (through a Ha emission criterion), and/or (2) through 
L-band or Spitzer 3.6 /im excess, and/or (3) their emission beyond 3.6 /im was classified 



1 Although originally quoted as 2.5-5 Myr, the timescale was recently updated by the same group ifliil 
using revised heat production values for 26 Al decay. 
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FIGURE 1. Age of stellar sample vs. fraction of stars with primordial disks (the "Haisch-Lada 2 " plot) 
either through Ha emission or infrared excess diagnostics. The best fit exponential decay curve is plotted 
with timescale T disk = 2.5 Myr. Disk fraction data are plotted for (in age order) NGC 2024 [0.3 Mvr;l29tl. 
NGC 1333 [1 Myr; El, Taurus [1.5 Myr; El, Orion Nebula Cluster [1.5 Myr; III, NGC 7129T2 Myr; 
IH, NGC 2068/71 [2 Myr; El, Cha I[2.6 Myr; El El, IC 348 [2.5 Myr; EE cr Ori [3 Myr; HI, NGC 
2264 132 Myr; El, Tr 37 [4.2 Myr; El, Ori OB lb [4 Myr; El, Upper Sco [5 Myrjpl, NGC 2362 [5 
MvnEl, y Vel [5 Myr; El, A Ori [5 Myr; El, r) Cha [6 Myr®, TW Hya J8 Myr; EE 25 Ori [8 Myr; 
ESEITnGC 7160 [11.8 Myr; El, /3 Pic [12 Myr;|H, UCL/LCC [16Myr;gl. 



by the authors as being likely due to primordial disk due to the SED shape or strength 
of the IR excess. The nature of some disks is unclear. Lada et al. 1I2II1 and others have 
identified stars with weak IR excesses whose nature as stars with either accretion disks 
of lower optical depth or simply warm dusty debris disks is at present ambiguous. Given 
the rarity of "transition disksj ' (disks with large inner holes), their inclusion or exclusion 
is usually within the disk fraction uncertainties, and will have negligible impact on this 
analysis. The fraction of stars in the transition phase has been noted to be very high 
in a young sample [e.g. ~1 Myr CrA; HfJ. I have not yet attempted to disentangle the 
effects of stellar mass in Fig.Q], so the reader should simply interpret the disk fractions as 
being most representative of the low-mass population of stars (<2 M ) as they dominate 
the stellar samples. I have omitted results for more distant clusters whose disk fraction 



A glossary for common terminology for young stellar objects with disks (the "diskionary") was recently 
compiled by Evans et al. [14311 . 



statistics were completely dominated by high mass stars (>l-2 M ). 

Ages and age uncertainties are taken from the published studies, however minimum 
age uncertainties of ±1 Myr or ±30% (whichever is greater) were adopted if uncer- 
tainties were not quoted. The usual caveats exist for the ages plotted in Fig. [Q There 
are significant differences in the ages estimated using different evolutionary tracks, and 
even as a function of mass for a single set of tracks [47], and none of the tracks have 
consistently matched predictions of masses with dynamically-constrained masses over 
the stellar mass spectrum. 

Fig. [T] demonstrates that any statements regarding the lifetimes of primordial disks 
need to be statistical in nature. Statements to the effect that "all disks disappear" by ~3 
or ~6 Myr are oversimplified assessments. I have decided to be provocative and plot 
an exponential function to fit the data. An exponential is convenient as it has a value 
of one at age zero (assuming that all stars are born with disks), trends towards zero as 
age increases (all primordial accretion disks eventually disappear), and one only needs 
to fit one parameter (the e-folding timescale tdisk)- When simultaneously minimizing 
the residuals in the disk fraction and age (Fig.Q]), I arrive at a best fit timescale of x^isk 
= 2.5 Myr. Unfortunately the best fit has reduced Xv — 2.5, suggesting that either (1) 
the uncertainties in the assessed disk fractions and ages are significantly underestimated 
(by factor of Vl^5 ~ 1.6), and/or (2) there are significant cluster-to-cluster differences 
in disk fraction decay time, and/or (3) an exponential function is simply inadequate for 
fitting the trend. As cluster-to-cluster variations have been demonstrated - especially 
at age ~5 Myr [2^, [25L |2vL 53], and cluster ages are especially uncertain, I suspect 



that these two factors are the primary causes of the high % 2 of the best fit. Removing 
individual clusters from the fit varies T^isk by <10%, which is probably a reasonable 
estimate of the precision in x^isk given our current knowledge of the ages of these star- 
forming regions. 

Recent Spitzer surveys have quoted disk fractions as a function of stellar mass in clus- 
ter samples. We can now look at these data in a different way, and combine the various 
results and quote a single metric (Tdisk) t0 more concisely summarize the observed trend 
and allow comparison between non-coeval samples. The exponential decay formalism is 
convenient for calculating e-folding times for various subsamples where we have sparse 
data available. If disk fraction evolves as an exponential decay, and a cluster of a given 
age is observed to have a disk fraction that is assumed to lie on that decay curve, one 
can estimate the primordial disk decay timescale from that subsample: 

^disk = -*/Mfdtik) (!) 

Where x is the age of the sample, and fn s k is the observed disk fraction. Using this 
technique, I will estimate preliminary values of x^isk f° r stellar subsamples segregated 
by stellar mass. 

What is the characteristic timescale for the primordial disks around brown dwarfs ? 
Results from early surveys identifying 3-4 jim excesses among small samples of substel- 
lar objects hinted that half of all disks were likely dispersed within ~l-3 Myr i49l l50ll. 
More recent disk fractions for larger samples of substellar objects have been quoted in 
a series of papers by Luhman and collaborators, notably for IC 348 [x = 2.5 Myr; f^ 



42%;|51D, Cha I [T ~ 2.5 Myr; f dlsk = 50%; 51], and o Ori [t = 3 Myr; f disk ~ 60% 



2711 . To this set, I include the disk fraction of substellar objects in Upper Sco (t = 5 Myr) 



from Mohanty et al. [50], which ranged from ~5-20% depending on whether the ac- 
cretors were defined via spectroscopic or photometric techniques (we assume 12.5 ± 7.5 
%). Note that the individual fd\ s k values have uncertainties of ~7-20%, which translate 
into significant errors in Tdisk, especially for small fduk- Using equation Q3 these disk 
fractions and ages translate into decay timescales of v^isk — 2.9, 2.9, 5.9, and 2.4 Myr, 
for IC 348, Cha I, o Ori, and Upper Sco, respectively. These results suggest that the typ- 
ical primordial disk decay timescale for brown dwarfs is approximately z^isk — 3 Myr - 
i.e. marginally longer than for stellar samples (~2.5 Myr). 

What is the characteristic decay timescale for primordial disks around high mass 
stars? The small disk fractions, and low numbers of high mass stars in stellar groups, 
make this surprisingly difficult to quantify. Hernandez et al. [|52ll conducted a systematic 
survey of the nearest OB associations to quantify the fraction of Herbig Ae/Be stars 
(>2M ). Their results found f^ < 5% for all of their samples (3-15 Myr), and for 
their two youngest samples: fduk = 5.1 ±2.0% in Ori OBlbc (3.5 ±3 Myr) and fdisk 
= 4.3 ± 1.8% in Tr 37. Using Spitzer to survey the ~2.5 Myr-old IC 348 group, Lada 
et al. ||2~Hl found a disk fraction of 1 1 ± 6% among > 1 .3 M stars. The Hernandez et al. 
0511 survey of the ~3 Myr-old o Ori clusters does not provide spectral types or masses, 
but interpolation of their Fig. 1 1 and Table 3 suggests a disk fraction of >1.3 M stars 
of ~10%, consistent with Lada et al.'s findings for IC 348. I122ll found no evidence for 
primordial disks around a sample of 92 BAFG-type (> 1 .3 M ) members of ~5 Myr-old 
Upper Sco (consistent with f^jt < 4%; 95% confidence), however 7/21 (~35%; K0-K6) 
of ~ 1-1.3 M stars have primordial disks! The results from these three Spitzer surveys 
are consistent with %disk — 1.2 Myr for >1.3 M stars. These results are also broadly 
consistent with the frequency of >2 M Herbig Ae/Be stars in nearby associations Q52Q . 

It is possible that tdisk could be used as a very coarse age estimator for multi- 
wavelength investigations of distant star-forming regions imaged both in the infrared 
and in X-rays where there is an estimate of the young stellar population both with 
primordial disks (class O/I/II objects) and without (class III objects). The mean age 
of the population would then be approximately x ~ -In fdisk x ^disk where fdisk = ^disk 
I (^disk + N„ disk) an d ^disk ~ 2.5 Myr. This estimate would only provide the coarsest 
of ages (as we now have evidence that disk lifetimes are dependent on stellar mass and 
environment), however in the quest for useful stellar age estimators, even the bluntest of 
diagnostics can be helpful. 

In summary, it is clear from Spitzer surveys that the lifetime of primordial disks is 
not only a function of age, but stellar mass, multiplicity, and proximity to O-type stars. 
Disk fraction appears to vary roughly as an exponential decay, with typical timescale 
Xdi s k — 2.5 Myr. This constant appears to vary from ~1.2 Myr for >1.3 M stars to ~3 
Myr for brown dwarfs. Although numerous mechanisms have been posited for depleting 
circumstellar disks, we need more observations to better constrain the disk evolution as a 
function of these stellar parameters (and for other untested parameters, e.g. metallicity), 
and more theoretical work to model these depletion mechanisms. It is clear that there 
are ample future opportunities for observations with the Subaru telescope to improve 
our understanding of the formation and early evolution of planetary systems. 
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